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Shunt  capacitors  are  commonly  used  in  distribution  system  for  reactive  power  compensation.  Different 
analytical,  numerical  programming,  heuristic  and  artificial  intelligent  based  techniques  have  been 
proposed  in  the  literature  for  optimum  shunt  capacitor  bank  (SCB)  placement.  This  paper  will  present 
a  very  detailed  overview  of  optimum  SCB  placement  techniques.  Six  different  approaches  of  optimum 
SCB  placement  based  on  minimization  of  power  losses,  weakest  voltage  bus  approach  and  maximization 
of  system  loadability  will  be  applied  on  four  different  radial  distribution  test  systems.  The  results  will  be 
compared  on  the  basis  of  power  loss  reduction,  voltage  profile  improvement,  system  loadability 
maximization  and  the  line  limit  constraint. 

©  2013  Elsevier  Ltd.  All  rights  reserved. 


Contents 


1.  Introduction . 429 

1.1.  Maximizing  the  distribution  system  efficiency . 430 

1.2.  Reactive  power  management  in  deregulated  power  market . 430 

1.3.  Reactive  power  management  in  presence  of  distributed  generation . 431 

2.  Optimum  shunt  capacitor  placement  techniques— A  review . 431 

2.1.  Analytical  methods . 431 

2.2.  Numerical  programming  methods . 432 

2.3.  Heuristics  methods . 432 

2.4.  Artificial  intelligent  methods . 432 

2.5.  Multi-dimensional  problems . 433 

3.  Comparitive  study . 433 

3.1.  Simulation  and  results . 434 

4.  Discussion . 434 

5.  Conclusion . 435 

Acknowledgment . 437 

Appendix . 437 

References . 437 


1.  Introduction 

Prior  to  1950s  the  shunt  capacitor  banks  (SCB)  were  placed 
nearer  to  the  main  substation  for  capacitive  reactive  power 
compensation,  it  helps  in  improving  the  power  factor,  reduces 


*  Corresponding  author.  Tel.:  +60  3  79675238;  fax:  +60  3  79675316. 
E-mail  address:  mohsinaman@gmail.com  (M.M.  Aman). 

1364-0321  /$  -  see  front  matter  ©  2013  Elsevier  Ltd.  All  rights  reserved. 
http://dx.doi.Org/10.1016/j.rser.2013.10.002 


I* 1 2R  power  losses  and  improving  the  voltage  profile.  SCB  changes 
the  power  losses  up  to  the  point  of  coupling,  however  to  get  the 
maximum  benefit  it  must  be  placed  as  nearer  to  the  load  as 
possible.  With  the  availability  of  pole  mounted  equipment  includ¬ 
ing  SCB,  the  trend  has  changed.  The  capacitor  banks  are  now 
placed  on  primary  distribution  lines  as  well  [1-5], 

The  capacitor  unit  is  considered  as  the  basic  building  block  of  SCB. 
Capacitor  units  are  connected  in  paralleled-series  combinations  and 
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form  a  single-phase  capacitor  bank,  within  a  steel  enclosure.  The 
series  combination  reduces  the  cost  of  dielectric  while  parallel 
combination  increase  the  total  capacitance  of  SCB.  As  a  general  rule, 
the  minimum  number  of  units  connected  in  parallel  is  such  that 
isolation  of  one  capacitor  unit  in  a  group  should  not  cause  a  voltage 
unbalance  more  than  110%  of  rated  voltage  on  the  remaining 
capacitors  of  the  group.  Equally,  the  minimum  number  of  series 
connected  groups  is  that  in  which  the  complete  bypass  of  the  group 
does  not  subject  the  others  remaining  in  service  to  a  permanent 
overvoltage  of  more  than  110%  [2].  The  amount  of  reactive  power 
(Qc)  from  capacitor  depends  on  applied  voltage  (V)  and  capacitive 
reactance  (Xc),  given  by  Eq.  (1)  [6], 

Q_c  =  V2/Xc  (1) 

The  recent  power  system  blackouts  [7,8]  due  to  insufficient 
reactive  power  have  also  resulted  in  focused  towards  meeting 
reactive  power  demand  of  the  system  locally  using  static  capacitor 
banks.  The  combined  US  Canada  task  force  on  August  2004 
blackout  also  concluded  that  the  reactive  power  supplies  in 
Northeast  Ohio  were  exhausted  which  resulted  in  loss  of  several 
critical  bulk  power  supply  systems  and  helped  cascaded  generator 
interruptions  [9],  During  low  voltage  emergencies  e.g.  generator 
rescheduling,  line  restoration  or  operated  directed  load  tripping, 
the  author  in  [10,6]  proposed  shunt  capacitor  bank  series  group 
shorting  (CAPS)  method.  CAPS  shorted  section  increases  the 
reactive  power  supplied  during  periods  of  low  voltages  by  shorting 
several  series  groups  of  capacitor  units  (Qct  =  V2/Xcl).  The 
shorted  section  in  CAPS  comprises  of  20%  to  33%  of  total  bank. 
The  detailed  study  and  feasibility  of  CAPS  on  EHV  and  HV  network 
are  also  presented  in  [11],  In  case  of  highly  loaded  systems,  it  is 
believed  that  the  optimum  capacitor  placement  solve  the  mini¬ 
mization  of  losses  more  adequately  and  optimum  setting  of 
voltage  regulators  solve  the  voltage  drop  problems  in  a  better 
manner  [3], 

The  need  for  reactive  power  support  in  distribution  system 
may  be  arises  due  to  the  following  reasons. 

3.3.  Maximizing  the  distribution  system  efficiency 

Distribution  system  usually  suffers  from  two  major  problems, 
high  power  losses  and  poor  voltage  profile.  Losses  are  defined  as 
the  difference  between  the  energy  into  the  system  and  the  energy 
that  is  utilized  by  the  end  users.  Generally  electric  system  losses 
can  be  categorized  as  technical  or  non-technical  losses  [12,13], 
Technical  losses  in  distribution  system  occurs  at  different  stages 
from  the  main  substation  till  the  consumer  end,  including  substa¬ 
tion  transformer,  primary  lines,  line  equipment  voltage  regulators 
and  surge  arrestor,  distribution  transformer,  secondary  lines  and 
consumer  services.  The  loss  calculation  methods  at  different  stages 
are  discussed  in  detail  in  [14],  Electric  Power  Research  Institute 
(EPRI)  and  Energy  Information  Administration  (EIA)  of  U.S.  con¬ 
cluded  that  [14,15]: 


1.  The  distribution  losses  range  from  33.7%  to  64.9%  of  the  total 
system  losses. 

2.  About  7%  of  the  total  electricity  production  is  transmitted  in  the 
United  States  as  transmission  and  distribution  losses  [15], 

3.  EPRI  research  on  42  distribution  circuits  estimated  that  54%  of 
total  losses  are  occurring  in  distribution  transformer  (although 
the  efficiency  of  distribution  transformer  lies  above  99%)  and 
38%  of  total  distribution  losses  are  occurring  in  primary  lines,  as 
shown  in  Fig.  1 .  One  of  the  major  reason  for  such  a  high  losses 
in  distribution  transformer  is  total  number  of  distribution 
transformer  placed  in  distribution  system  [14],  In  2003,  it  is 
estimated  that  there  are  50  million  distribution  transformers  in 
use  in  United  States  [16], 

With  such  a  higher  power  losses  in  distribution  system,  it  is 
highly  necessary  to  reduce  the  line  losses  occurring  in  primary  line 
as  much  as  possible.  The  higher  losses  results  in  limiting  the  line 
capacity  ( thermal  limits )  as  well  as  higher  voltage  drop  ( voltage 
limits )  in  the  power  system.  In  literature  it  has  also  been  concluded 
that  the  maximum  loading  of  the  distribution  system  is  limited  by 
the  voltage  limit  rather  than  the  thermal  limit  [17].  Large  power 
consumers  also  installed  shunt  capacitor  to  improve  the  overall 
power  factor  and  thus  save  the  cost  of  poor  power  factor  penalty. 
Three  different  compensation  techniques  are  available  in  literature 
including  individual  compensation;  group  compensation  and 
centralized  compensation  to  improve  the  power  factor.  Any  one 
of  the  method  or  all  of  the  method  can  be  utilized  to  take  the 
maximum  advantage  of  improved  power  factor  [18].  For  power 
factor  improvement,  an  unloaded  synchronous  motor  can  also  be 
used  instead  of  shunt  capacitor.  The  amount  of  reactive  power  is 
controlled  from  its  excitation  system  and  thus  it  behaves  like  a 
variable  capacitor  [19].  Now-a-days,  manufacturers  are  bound  to 
design  electrical  equipment  with  higher  efficiency  and  high  power 
factor  [20], 

3.2.  Reactive  power  management  in  deregulated  power  market 

With  the  restructuring  of  power  system,  the  complexity  of 
power  system  has  been  increased.  The  vertically  integrated  power 
system  has  been  separated  into  GenCOs,  TransCOs,  and  DisCOs. 
A  central  regulating  company  Independent  System  Operator  (ISO) 
and  Regional  Transmission  Organizations  (RTO)  has  been  formed 
for  maintaining  the  quality,  reliability  and  security  of  electric 
service  [21,22],  It  is  also  a  fact  that  the  existing  transmission 
systems  in  most  of  the  countries  are  quite  old.  For  example  in 
United  States,  the  345  kV  bulk  transmission  system  and  associated 
substation,  cables  and  wires  are  40  years  old  and  above  [23].  Such 
a  system  is  not  able  to  meet  the  growing  demand  and  transfer  the 
generated  power  from  the  centralized  generation  to  the  distribu¬ 
tion  system.  Transmission  investment  has  been  falling  for  a 
quarter  century  at  an  average  rate  of  almost  US$50  million  a  year 
(in  constant  2003  U.S.  dollars),  however  there  has  been  a  small 
upturn  in  the  last  few  years  [23].  Other  than  constructing  new 
transmission  lines  there  are  other  options  to  release  the  transmis¬ 
sion  system  congestion  including  distributed  generator  placement, 
static  capacitor  bank,  FACTS  (Flexible  AC  Transmission  Systems) 
devices,  voltage  regulators  and  energy  conservation  [24,25],  It  is 
also  a  fact  that  the  generation  the  addition  of  extra  kVAR  on  a 
generator  operating  at  0.9  power  factor  decreases  the  amount  of 
real  power  output  by  about  half  a  kilowatt  [26],  Thus  in  restruc¬ 
tured  power  market,  the  generator  companies  prefer  to  generate 
maximum  active  power  (kW)  and  get  maximum  profit  ($). 

The  reactive  power  (kVAR)  market  is  not  as  simple  as  real 
power  (kW)  market.  The  fundamental  difficulty  with  reactive 
power  markets  is  that  reactive  power  does  not  “travel”  far,  thus 
it  is  expected  there  would  be  extreme  local  geographical  market 
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power  in  the  provision  of  reactive  power  [27].  This  philosophy  of 
reactive  power  market  suggests  the  following  things  [26]: 

1.  The  generator  should  maintain  close  to  unity  power  factor 
under  normal  conditions  using  capacitor  bank.  However  gen¬ 
erator  dynamic  response  is  used  to  provide  reactive  reserves  in 
case  of  contingency. 

2.  The  distribution  system  must  have  enough  reactive  power 
support  so  that  distribution  feeders  have  a  net  unity  power 
factor  or  even  slightly  leading  power  factor  at  low  loads  and  so 
that  transmission  system  steady  state  reactive  power  is  essen¬ 
tially  all  provided  by  transmission  system  capacitors. 


1.3.  Reactive  power  management  in  presence  of  distributed 
generation 

Electric  power  generation  that  is  integrated  within  the  dis¬ 
tribution  systems  are  known  as  “distributed”  or  “dispersed” 
generation  (DG).  The  DG  source  or  technology  can  be  a  traditional 
combustion  generator  (such  as  diesel  reciprocating  generator  and 
natural  gas-turbine)  and  non-traditional  generator  including  fuel 
cell,  storage  device  and  renewable  energy  source  (such  as  wind 
turbine  and  photovoltaic).  Non-utility  companies  are  investing  in 
distribution  system  to  meet  the  active  power  demand  (MW)  and 
get  the  maximum  profit.  For  example  in  United  States  the 
percentage  of  nonutility  generators  (NUGs)  has  increased  from 
40GW  to  more  than  150  GW  in  10  years  from  1990  to  2000 
[28-32], 

DCs  are  mainly  considered  as  active  source  of  energy  [33], 
however  at  higher  system  loading  with  maximum  DG  penetration, 
the  poor  voltage  profile  can  be  a  big  challenge  for  the  system 
operator  thus  the  reactive  power  compensation  approach  must  be 
utilized  to  maintain  the  voltages  in  allowable  limits  [34], 
In  restructured  power  system,  the  system  operator  is  responsible 
to  maintain  the  allowable  voltages  by  call  upon  the  GENCOs  to 
produce  more  or  less  reactive  power,  by  adjusting  the  field  current, 
by  adding  or  removing  reactive  power  devices  (capacitor  and 
reactors)  or  changing  the  tap  changer  position  on  the  transformer. 
The  author  in  [35]  has  analyzed  the  importance  of  reactive  power 
in  presence  of  DG  and  conclude  that  the  presence  of  DG  may 
results  in  voltage  rise  problem  at  light  load,  thus  the  voltage 
regulating  device  must  be  also  presented.  The  energy  curtailment 
from  DG  is  not  a  good  solution  as  this  will  result  in  revenue  lost. 
DG  presence  in  the  system  also  affects  the  reactive  power  manage¬ 
ment  plan.  For  example  in  case  of  wind  generation,  asynchronous 
induction  generators  are  used.  Such  generators  need  reactive 
power  from  the  system  to  which  they  are  connected.  Different 
methods  of  reactive  compensations  are  stated  in  literature  [36] 
including  synchronous  generator,  shunt  capacitor  banks  and  end- 
user  reactive  power  compensation  within  the  reactive  power 
consumption  equipment.  Here  it  is  also  need  to  be  mentioned 
that  the  growing  trend  of  using  non-conventional  power  genera¬ 
tion  (using  wind  and  solar  energy)  has  led  to  the  bounding  that 
the  renewable  energy  generation  must  also  play  their  role  in 
improving  the  voltage  profile  and  providing  necessary  reactive 
power  support.  Now-a-days  state  of  the  art  technology  has  come 
out,  the  wind  generation  is  now  using  doubly  fed  induction 
generator  and  PV  inverters  are  using  special  self-commutated  line 
inverter,  capable  of  absorbing  and  supplying  reactive  power  at 
different  system  loading.  The  reactive  power  capability  of  solar 
and  wind  power  plants  can  be  further  enhanced  by  the  addition  of 
SVC,  STATCOMS  and  other  reactive  support  equipment  at  the  plant 
level.  Currently,  inverter-based  reactive  capability  is  more  costly 
compared  to  the  same  capability  supplied  by  synchronous 
machines  [37-39], 


Thus  the  importance  of  reactive  power  compensation  has  been 
significantly  increased  due  to  restructuring  of  power  market  and 
conversion  of  passive  network  to  active  network.  Utilities  are  also 
focussing  towards  reduction  in  power  losses  and  utilizing  the 
distribution  lines  efficiently.  SCB  being  the  low  initial  cost,  no 
maintenance  and  no  personnel  cost  is  the  most  cost  effective 
solution  for  reactive  power  compensation.  Section  2  will  present  a 
detailed  overview  of  capacitor  bank  placement  techniques  stated 
in  the  literature.  The  non-optimum  placement  or  sizing  of  SCB 
may  result  in  increased  power  losses  as  losses  versus  capacitive 
MVAr  follows  the  deep  bath  curve  relationship  [40].  Thus  it  is 
necessary  to  optimally  place  VAR  equipment  in  the  distribution 
system. 

2.  Optimum  shunt  capacitor  placement  techniques— A  review 

In  literature,  different  authors  have  proposed  different  methods 
considering  different  fitness  function  including  minimization  of 
power  losses,  reduction  in  installation  cost,  improvement  in 
voltage  profile,  lessen  the  burden  on  existing  lines,  maximization 
of  system  stability  and  others.  Shunt  capacitors  are  placed  in  a 
combination  of  fixed  and  switched  (variable)  capacitor  banks.  The 
size  of  fixed  capacitor  bank  depends  on  average  reactive  power 
demand  of  the  system  while  switched  capacitor  banks  supply 
difference  of  current  reactive  power  demand  and  fixed  capacitive 
power  available.  Special  control  mechanism  is  used  to  control  the 
switched  capacitor  bank  power.  Shunt  capacitors  (Qc)  are  available 
in  market  in  discrete  sizes  which  are  multiple  (k)  of  smallest 
capacitor  size  (Q-nm),  given  by  Eq.  (2)  [41,42],  However  authors 
have  proposed  such  methods  which  results  in  both  continuous 
and  discrete  capacitor  size.  In  case  of  continuous  capacitor  size,  it 
will  be  assumed  that  the  capacitor  size  will  be  a  combination  of 
fixed  and  switched  capacitor.  The  coarse  tuning  of  capacitor  under 
steady  state  condition  will  be  achieved  using  switched  capacitor 
bank. 

Qc  =  k  X  Qmin  (2) 

The  capacitor  placement  problems  from  the  literature  can  be 
categorized  into  analytical  methods,  numerical  programming 
methods,  heuristic  methods  and  artificial  intelligent  methods. 
Authors  have  also  combined  the  capacitor  problem  with  other 
power  system  problems  including  distributed  generation,  network 
reconfiguration  and  voltage  regulator. 

2.1.  Analytical  methods 

In  earlier  time,  when  powerful  computing  resources  were 
unavailable  or  expensive,  the  author  proposed  calculus  based 
analytical  algorithm.  The  authors  have  also  made  some  approx¬ 
imation  in  order  to  reduce  the  computation  procedure.  Analytical 
methods  have  also  been  proposed  in  optimum  capacitor  place¬ 
ment  and  sizing.  The  initial  work  was  carried  out  by  Neagle  [1]  in 
1956  for  optimum  single  and  multiple  capacitor  bank  in  case  of 
uniform  and  non-uniform  distribution  of  load.  He  suggested  that 
in  uniformly  distributed  load  the  capacitor  bank  must  be  placed  at 
1  -(1/2)  (capacitive  kvar/system  kvar)  distance  from  the  main 
substation.  Cook  worked  on  the  same  guideline  of  [1]  and 
proposed  a  more  practical  algorithm  for  fixed  capacitor  bank  for 
uniformly  distributed  considering  average  reactive  load  in  the 
system  [43].  Cook  [43]  suggested  that  the  optimum  location  of 
SCB  must  be  2/3  (Reactive  load  factor).  Cook  extended  his  work  in 
[44]  to  include  switched  capacitors.  Later  on  several  other  analytical 
based  methods  have  been  proposed  for  the  capacitor  placement 
[45,46],  Schmill  [47]  extended  the  work  of  Cook  [43],  Equations  are 


432 


M.M.  Aman  et  al.  /  Renewable  and  Sustainable  Energy  Reviews  30  (2014)  429-439 


given  for  sizing  and  placement  of  n  capacitors  on  a  uniform  feeder 
with  a  uniformly  distributed  load.  The  necessary  conditions  for 
optimal  sizing  and  placement  of  one  or  two  capacitors  on  a  feeder 
with  discrete  loads  and  non-uniform  resistance  are  presented.  An 
iterative  approach  is  suggested  to  solve  the  problem.  Chang  et  al. 
[48,49]  assumes  a  feeder  with  a  uniform  load  and  a  concentrated 
end  load.  Accounting  for  both  peak  power  losses  and  energy 
losses,  he  determines  the  optimal  location  of  a  fixed  capacitor 
and  the  resulting  savings,  given  the  capacitor  size.  The  optimal 
solution  is  determined  by  considering  each  of  the  available 
capacitor  sizes. 

Analytical  methods  are  considered  as  simple  methods,  however 
implication  of  such  methods  to  solve  the  capacitor  placement 
problem  required  all  assumptions  (e.g.  load  variation  neglected) 
and  scenarios  (e.g.  distributed  non-distributed)  as  considered  by 
the  author  in  developing  that  algorithm.  The  author  in  [50]  has 
indicated  that  the  famous  “2/3  rule”  may  results  in  negative  saving 
if  considered  in  different  scenario.  Here  it  is  also  noted  that  most 
of  the  analytical  methods  discussed  earlier  considered  modelling 
of  the  capacitor  placement  locations  and  sizes  as  continuous 
variables.  Therefore,  the  results  would  need  to  be  rounded  up  or 
down  to  the  nearest  practical  value,  which  may  result  in  an 
overvoltage  problem  or  loss  savings  ($)  less  than  the  calculated 
one.  The  more  recent  analytical  methods  [46,51-54]  are  much 
more  accurate  and  practical  for  distribution  systems. 

2.2.  Numerical  programming  methods 

Numerical  Programming  is  a  technique  by  which  mathematical 
problems  are  formulated  so  that  they  can  be  solved  with  arith¬ 
metic  operations.  Numerical  programming  methods  are  iterative 
techniques  used  to  maximize  (or  minimize)  an  objective  function 
of  decision  variables.  The  values  of  the  decision  variables  must  also 
satisfy  a  set  of  constraints.  With  the  availability  of  fast  computing 
skills  and  large  memory  availability,  the  utilization  of  numerical 
methods  in  power  system  has  increased  [55,56],  In  optimum 
capacitor  location  problem,  authors  have  formulated  different 
mathematical  models  and  have  utilized  numerical  methods  to 
solve  the  problem.  Duran  in  [57],  have  used  the  dynamic  pro¬ 
gramming  approach  and  implemented  the  Schmil  work  [47]  of 
uniformly  and  randomly  distributed  load  to  find  the  optimum 
capacitor  placement.  The  formulation  in  [57]  is  simple  and  only 
considers  the  energy  loss  reduction  and  accounts  for  discrete 
capacitor  sizes.  Fawzi  et  al.  [58]  extended  the  work  of  [57]  and 
included  the  released  kV  A  into  the  savings  function.  Ponnavsikko 
and  Rao  [59]  used  a  numerical  method  called  the  method  of  local 
variations  and  further  expanded  the  problem  to  include  the  effects 
of  load  growth,  and  switched  capacitors  for  varying  load.  Lee  [50] 
has  proposed  iterative  based  optimization  technique  considering 
net  monetary  saving  as  convergence  criteria  in  optimum  place¬ 
ment  of  fixed  and  switched  capacitor  banks.  Later  on  Baran  and 
Wu  [60,61]  used  the  mixed  integer  programing  approach  to  solve 
the  capacitor  problem.  Sharaf  used  the  full  load  flow  model  to  find 
the  optimum  shunt  capacitor  placement  in  distribution  system.  He 
mentioned  that  the  equivalent  model  proposed  in  [62]  is  not 
useful  in  finding  optimum  SCB  location  as  the  receiving-end 
voltage  on  a  distribution  system  decreases  quadratically  as  system 
load  increases  [63].  Khodr  [64]  solved  the  SCB  problem  using 
mixed-integer  linear  problem  and  considered  overall  energy  sav¬ 
ing  as  fitness  function. 

2.3.  Heuristics  methods 

Heuristics  methods  are  “hints”,  “suggestions”,  or  “rules  of 
thumb”  based  methods  that  are  developed  through  intuition, 
experience,  and  judgment.  Heuristic  methods  produce  fast  and 


practical  strategies  which  reduce  the  exhaustive  search  space  and 
can  lead  to  a  solution  that  is  nearer  to  the  optimal  solution  with 
confidence  [65,66],  Heuristic  based  techniques  are  also  commonly 
used  to  solve  shunt  capacitor  placement  problem  [3,67-71], 
Abdel-Salam  et  al.  [67]  proposed  a  heuristic  technique  based  on 
identifying  the  sensitive  node  and  placing  such  capacitor  bank 
(size)  giving  the  greatest  loss  reduction  due  to  capacitor  place¬ 
ment.  Chis  et  al.  in  [68]  extended  the  work  of  [67]  and  considered 
cost  of  capacitor  bank  (size)  from  both  energy  and  peak  power  loss 
reductions.  In  [69],  the  author  has  evaluated  bus-bar  sensitivity 
index  to  decide  the  capacitor  position(s),  sigmoid  function  is 
utilized  to  find  the  discrete  value  of  capacitor  size.  In  latest 
research  [70],  the  author  has  used  node  voltage  stability  index  to 
find  the  candidate  bus  for  capacitor  position  and  maximization  of 
net  savings  from  power  loss  reduction  and  the  capacitors  invest¬ 
ment  as  optimum  capacitor  size.  In  [40],  the  author  has  considered 
the  end  of  the  weakest  line  as  the  potential  candidate  for  optimum 
shunt  capacitor  placement  and  capacitor  size  is  selected  corre¬ 
sponding  to  the  minimum  power  losses  using  PSO  algorithm. 
In  [71],  the  author  has  used  direct  search  algorithm  to  determine  the 
optimal  sizes  of  fixed  and  switched  capacitors  together  with  their 
optimal  locations  in  a  radial  distribution  system  so  that  net  savings 
are  maximized  and  improvement  in  the  voltage  profile  is  achieved. 

2.4.  Artificial  intelligent  methods 

The  simplest  of  search  algorithms  in  optimization  is  exhaustive 
search  that  tries  all  possible  solutions  from  a  predetermined  set 
and  subsequently  picks  the  best  one.  However  such  methods  are 
considered  as  non-efficient  in  terms  of  computation  time  and 
space  requirement.  Intelligent,  greedy  and  nature  observed  heur¬ 
istic  techniques  have  also  been  proposed  in  literature  [72], 
commonly  known  as  artificial  intelligent  (Al)  methods.  Al  methods 
are  a  special  class  of  heuristic  search  methods  [72].  Intelligent 
based  optimization  methods  have  also  been  utilized  in  finding 
optimum  shunt  capacitor  location  and  sizing.  Authors  have  used 
genetic  algorithms  (GA)  [3,73-77]  ,  fuzzy  [78,79],  immune  algo¬ 
rithm  [80,81],  Tabu  search  [82],  fuzzy-GA  [83],  particle  swarm 
optimization  [84],  plant  growth  simulation  algorithm  [41,85], 
memetic-algorithm  approach  [86],  teaching  learning  based  opti¬ 
mization  algorithm  [87],  ant  colony  [88],  graph  search  algorithm 
[89]  and  hybrid  algorithm  [90-92], 

In  [93],  the  author  has  used  non-dominated  sorting  genetic 
algorithm  (NSGA)  to  solve  multi-criteria  problem  in  capacitor 
placement.  A  multi-criteria  optimization  problem  requires  simul¬ 
taneous  optimization  of  a  number  of  objectives  with  different 
individual  optima  [11,14],  Objectives  are  such  that  none  of  them 
can  be  improved  without  degradation  of  another.  Hence,  instead  of 
a  unique  optimal  solution,  there  exists  a  set  of  optimal  trade-offs 
between  the  objectives,  the  so-called  pareto-optimal  solutions. 

The  author  in  [94,95]  has  presented  an  algorithm  for  optimiz¬ 
ing  shunt  capacitor  sizes  on  radial  distribution  lines  with  non- 
sinusoidal  substation  voltages,  such  that  the  rms  voltages  and 
their  corresponding  total  harmonic  distortion  lie  within  prescribed 
values.  The  result  shows  that  the  optimal  capacitor  sizes  found  by 
neglecting  the  harmonic  components  may  result  in  unacceptable 
voltage  distortion  levels.  In  [96],  the  author  has  considered  the 
presence  of  non-linear  load  in  distribution  system  in  solving  the 
optimal  capacitor  placement  problem.  The  authors  has  used  PSO  to 
search  for  optimal  locations,  types,  and  sizes  of  capacitors  to  be 
placed  and  optimal  numbers  of  switched  capacitor  banks  at 
different  load  levels.  In  [97],  the  author  has  solved  the  problem 
of  capacitor  placement  using  genetic  algorithm  (GA)  in  the  pre¬ 
sence  of  voltage  and  current  harmonics.  The  author  has  claimed 
that  the  proposed  method  results  in  lower  THDV  and  greater 
annual  benefits  as  compared  to  [98-100],  Presently  Al  methods  are 
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Table  1 

Criterion  for  comparing  the  different  methods  of  SCB  Placement. 


Parameters 

Formulae 

Abbreviations 

Total  Line  Loss  Reduction  (TLLR) 

Voltage  Profile  Improvement  (VPI) 

System  Maximum  Loading  Improvement 
(SMLI) 

Line  Limit 

Re{lo5se5}0-Re{losses}se, 

Re  { losses  }0 

ypjo,  ..  1 0O 

vu/°- - os-.va),, - x  1UU 

SML1%=  x  ioo 

4max(0) 

■  Subscript  (0)  is  representing  the  base  case  when  no  SCB  is  present  in  the  system. 

■  Subscript  (SCB)  is  representing  the  case  when  SCB  is  present  in  the  system. 

■  Amax  is  maximum  loadability  of  the  system. 

■  Vi  is  voltage  magnitude  at  bus  i. 

■  Li  is  active  load  at  bus  i  pu. 

lSkCB>I0k 

where  l<k<nbr 

■  nbus  is  total  number  of  buses. 

■  Ik  is  branch  current  in  kt h  branch. 

■  nbr  is  total  number  of  branches. 

considered  as  the  most  powerful  methods  in  solving  many  power 
system  problems  including  single-objective-multi-constraints,  multi- 
objective-multi-constraints  or  multi-criteria-pareto-optimal  pro¬ 
blems.  However  in  large  system,  AI  methods  suffers  from  high 
computation  time  and  large  memory  space  requirement,  in  such 
cases  hybrid  methods  are  considered  more  powerful. 

2.5.  Multi-dimensional  problems 

In  multi-dimensional  problems,  the  authors  have  combined  the 
capacitor  placement  problems  with  other  power  system  problems 
including  network  reconfiguration  [101-107],  DG  placement  [108- 
111],  voltage  regulators  [112-118]  and  load  tap  changer  [119],  The 
author  in  [112-114]  considered  the  capacitor  placement  problem  a 
multi-dimensional  problem  and  incorporated  the  voltage  regula¬ 
tor  placement  in  addition  to  control  the  volt/var.  The  author  in 
[120]  has  combined  different  type  of  DCs,  shunt  capacitor  and 
network  reconfiguration  in  a  single  problem  and  stated  that  power 
losses  are  significantly  reduced  when  all  three  objectives  are 
solved  simultaneously. 

In  next  section  (Section  3),  six  different  methods  based  on 
minimization  of  power  losses,  weakest  voltage  bus  and  maximiza¬ 
tion  of  system  loadability  for  optimum  shunt  capacitor  placement 
will  be  presented  and  discussed  in  detail. 

3.  Comparitive  study 


Method  2:  Grid  search  algorithm  [122] 

In  grid  search  algorithm,  shunt  capacitor  is  added  to  each  bus 
and  the  size  of  capacitor  is  changed  from  0%  to  100%  of  total 
connected  load  in  small  steps.  Minimization  of  losses  is 
considered  as  the  objective  function.  For  this  purpose,  succes¬ 
sive  load  flow  methods  are  used  for  each  step  of  capacitor  size. 
The  main  constraints  are  to  restrain  the  maximum  capacitor 
size  selected  as  total  load  size. 

Method  3:  Golden  section  search  algorithm  [122] 

In  golden  section  search  algorithm,  the  searching  space  is 
decreased  by  checking  some  discrete  values  of  capacitor  size 
only,  in  every  iteration.  Minimization  of  losses  is  considered  as 
an  objective  function.  The  main  constraints  are  to  restrain  the 
maximum  capacitor  size  selected  as  total  load  size. 

Methods  1-3  is  implemented  using  the  MATLAB  based  voltage 
stability  and  optimization  programming  (VS&OP)  tool  [122], 

Method  4:  Optimum  shunt  capacitor  placement  and  sizing  based 
on  min  of  losses  using  PSO 

In  this  method,  the  optimum  location  and  optimum  size  is 
determined  simultaneously  using  particle  swarm  optimization 
(PSO)  algorithm.  PSO  has  advantage  over  other  methods  in  terms 
of  simplicity,  fast  convergence  as  compared  to  other  algorithms  (e. 
g.  genetic  algorithm)  [123,124],  The  brief  overview  of  PSO  method 
is  presented  in  the  appendix.  Minimization  of  active  PR  power 
losses  (PL),  given  by  Eq.  (5)  is  considered  as  fitness  function. 


In  this  section,  six  different  methods  for  optimum  shunt 
capacitor  placement  and  sizing  are  compared  on  the  basis  of  power 
losses  reduction,  voltage  profile  improvement,  maximization  of 
system  loadability  and  line  limit  constraint,  summarized  in  Table  1. 

Following  constraint,  given  by  Eqs.  (3)  and  (4)  is  used  through¬ 
out  the  analysis: 

Position  of  SCB:  2  <  QSCB  Position  <  nbus  (3) 

nbus 

Size  of  SCB:  QSCB  <  £  Q,  (4) 

i  =  2 

where  nbus  is  total  number  of  buses  and  Qi  is  the  reactive  power 
load  at  bus  i. 

The  brief  overview  of  methods  is  discussed  below: 

Method  1:  Analytical  method  [121] 

In  analytical  approach,  a  loss  sensitivity  factor  (dPtoss/dP,-)  is 
formulated  for  the  determination  of  the  optimum  size  and 
location  of  capacitor  bank  to  minimize  total  power  losses.  Like 
other  methods,  analytical  method  does  not  use  admittance 
matrix,  inverse  of  admittance  matrix  or  Jacobian  matrix.  This 
method  is  based  on  the  equivalent  current  injection  that  uses 
the  bus-injection  to  branch-current  (BIBC)  and  branch-current 
to  bus-voltage  (BCBV)  matrices. 


/  =  Min 


nbr 

I  W  Rf 


(5) 


l  i  =  1  J 

where  I  is  the  line  current,  R  is  the  resistance  of  branch  i.  nbr  is  total 
number  of  branches. 


Method  5:  Hybrid  method  based  on  weakest  line  and  minimiza¬ 
tion  of  power  system  losses  using  particle  swarm  optimization 
In  this  method,  the  SI  index  proposed  in  [125]  is  utilized  to  find  the 
weakest  voltage  bus  in  power  system  which  could  lead  to  voltage 
instability,  when  the  load  will  cross  the  critical  limit.  The  value  of 
index  is  given  by  Eq.  (6)  and  termed  as  stability  index  (SI). 

Sl(r)  =  \Vs\4—4  x  {Prx,j  — Qsr,j}2  — 4  x  {Prr,j  +  Qrx0)2  x  |VS|2  >  0  (6) 

where  l/s  is  sending  end  voltage,  Pr  is  active  load  at  receiving  end, 
Qr  is  reactive  load  at  receiving  end,  is  resistance  of  the  line  i-j,  x 
(ij)  is  reactance  of  the  line  i-j. 

Under  stable  operation,  the  value  of  SI  should  be  greater  than  zero 
for  all  buses.  When  the  value  of  SI  becomes  closer  to  one,  all  buses 
become  more  stable.  In  this  algorithm,  SI  value  is  calculated  for 
each  bus  in  the  network  and  sort  from  highest  to  lowest  value.  For 
the  bus  having  the  lowest  value  of  SI,  the  SCB  will  be  placed  at  that 
bus.  Once  the  SCB  location  has  been  identified,  the  size  of 
capacitor  bank  is  calculated  based  on  minimization  of  active  PR 
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Table  2 

Base  case  (before  capacitor  placement). 


Test 

system 

P-loss 

(MW) 

System 

loadability 

Vmj 

'^min 

Sum 

Mb) 

Weakest  voltage 
bus 

12-Bus 

0.0207 

5.32 

1/0.9439 

0.4204 

12 

30-Bus 

0.8819 

2.79 

1/0.8825 

8.0290 

27 

33-Bus 

0.2110 

3.41 

1/0.9038 

3.5228 

18 

69-Bus 

0.2250 

3.22 

1/0.9092 

3.6191 

65 

Table  3 

Capacitor  placement  based  on  different  methods. 

Method 

Capacitor  position/ 
size  (MVAr) 

P-loss 

(MW) 

System 

loadability 

V,„.«/ 

l^min 

Sum  (ViU) 

Method  1: 

Analytical  method 

12-Bus 

9/0.2103 

0.0126 

5.65 

1/0.9563 

0.4236 

30-Bus 

22/3.13572 

0.6831 

2.94 

1/0.8979 

8.1141 

33-Bus 

30/1.2298 

0.1514 

3.6 

1/0.9162 

3.5651 

69-Bus 

61/1.2920 

0.1521 

3.49 

1/0.9302 

3.6582 

Method  2: 

Grid  search  algorithm 

12-Bus 

9/0.2106 

0.0126 

5.65 

1/0.9563 

0.4236 

30-Bus 

21/3.45818 

0.6812 

2.93 

1/0.8979 

8.1183 

33-Bus 

30/1.265 

0.1514 

3.6 

1/0.9165 

3.5662 

69-Bus 

61/1.3203 

0.1520 

3.49 

1/0.9306 

3.6590 

Method  3: 

Golden  section  search 

algorithm 

12-Bus 

9/0.2102 

0.0126 

5.65 

1/0.9563 

0.4236 

30-Bus 

21/3.4347 

0.6812 

2.93 

1/0.8978 

8.1178 

33-Bus 

30/1.258 

0.1514 

3.6 

1/0.9165 

3.5660 

69-Bus 

61/1.330 

0.1520 

3.49 

1/0.9307 

3.6592 

Method  4: 

Minimization  of  power  losses 

12-Bus 

9/0.2102 

0.0126 

5.65 

1/0.9563 

0.4236 

30-Bus 

21/3.4347 

0.6812 

2.93 

1/0.8978 

8.1178 

33-Bus 

30/1.2580 

0.1514 

3.6 

1/0.9165 

3.5660 

69-Bus 

61/1.330 

0.1520 

3.49 

1/0.9307 

3.6592 

Method  5: 

Hybrid  (voltage  stability + power  losses) 

12-Bus 

12/0.1737 

0.0134 

5.66 

1/0.9561 

0.4232 

30-Bus 

27/2.384 

0.7148 

2.95 

1/0.8971 

8.0956 

33-Bus 

18/0.4798 

0.1891 

3.56 

1/0.9211 

3.5471 

69-Bus 

65/0.9822 

0.1691 

3.42 

1/0.9280 

3.6496 

Method  6: 

Maximization  of  system  loadability 

12-Bus 

11/0.4300 

0.0263 

6.06 

1/0.9710 

0.4266 

30-Bus 

24/5.15 

0.8647 

3.06 

1/0.9086 

8.1562 

33-Bus 

12/2.260 

0.2733 

3.81 

1/0.9357 

3.6069 

69-Bus 

62/2.64 

0.2214 

3.72 

1/0.9493 

3.6923 

power  losses  (PL),  given  by  Eq.  (5)  using  PSO  optimization 
algorithm.  The  similar  approach,  based  on  weaest  voltage  bus 
has  been  utilized  in  [40]  and  [70], 

Method  6:  Optimum  shunt  capacitor  placement  and  sizing  based  on 
maximization  of  system  loadability  using  PSO  (proposed  method) 

In  this  method,  PSO  method  will  be  utilized  to  find  the  optimum 
SCB  position  and  optimum  size,  based  on  maximization  of  system 
loadability  ( Amax )  as  given  by  Eq.  (7). 

f  =  MaX{y Imaxl  (7) 

To  find  the  maximum  loadability  of  the  system  (2max),  the  active 
and  reactive  load  is  increased  on  all  buses  (using  Eqs.  8a  and  8b) 
with  equal  loading  factor  of  0.05,  till  the  divergence  is  observed  in 
load  flow  analysis.  The  author  in  [126]  has  used  this  method  for 
optimum  network  reconfiguration  in  distribution  system. 


P o  x  loading  factor(^) 

(8a) 

:  Qo  x  loading  factor(A) 

(8b) 

a 

1.04 

1.02 

's' 

3  l 
■a 

*g  0-98 

D£ 

a 

S  0.96 
§&  0.95 

i  0.94 

> 

0.92 

0.9 

123456789  10  11  12 

Bus  No. 

b 

1.08 

1.06 

^  1.04 

3 
a. 

“  1.02 

<u 


'I  0-98 
I  0-96 

S)  0.94 

?  0.92 

> 

0.9 
0.88 

10  20  30  40  50  60  69 
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Fig.  2.  Comparison  of  voltage  profile  using  different  methods  of  shunt  capacitor 
placement,  (a)  12-Bus  Test  System,  and  (b)  69-Bus  Test  System. 

where  A  is  a  loading  factor.  Pa  and  Qo  is  initial  active  and  reactive 
power  load,  connected  with  ith  bus.  Pnew  and  Qnew  is  final  active 
and  reactive  power  load,  connected  with  ith  bus. 

3.t.  Simulation  and  results 

The  proposed  method  is  tested  on  12-bus  [127],  30-bus  [128], 
33-bus  [129]  and  69-bus  [60]  radial  distribution  test  systems  in  all 
three  scenarios.  Thukaram  load  flow  method  [130]  is  used  to  carry 
out  the  power  flow  analysis  for  the  radial  distribution  system.  In 
base  case,  when  no  capacitor  is  place  in  the  system,  following 
results  have  been  obtained,  given  in  Table  2. 

Table  3  is  showing  the  results,  when  different  methods  of 
optimum  capacitor  placement  and  sizing  are  applied  on  test 
systems. 


4.  Discussion 

The  above  methods  are  compared  on  the  basis  of  different 
criterion  given  in  Table  1.  Here  it  is  interesting  to  note  that  the 
method  1  to  method  4,  approximately  give  the  same  results.  Thus 
these  methods  will  be  discussed  simultaneously  as  “loss  reduction 
methods".  Following  important  points  are  concluded: 

1.  There  are  significant  improvements  in  reduction  in  I2R 
power  losses,  voltage  profile  improvement  and  in  maximiza¬ 
tion  of  system  loadability  (Amax)  in  comparison  with  base 
case  (no  capacitor).  However  the  method  based  on  loadability 
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Fig.  3.  Comparison  of  six  different  methods  of  optimum  shunt  capacitor  place¬ 
ment.  (a)  Total  Line  Loss  reduction,  (b)  Voltage  Profile  Improvement,  and  (c)  System 
Maximum  Loadability  Improvement. 


Table  4 

Maximum  loadability  improvement  due  to  single  capacitor  placement  considering 
voltage  limits. 


Methods 

Capacitor 

position 

Capacitor 

size 

(MVAR) 

Maximum 

loading 

margin 

Maximum 
load  (kV  A) 

Base  case 

_ 

_ 

0.89 

529.91 

Analytical 

9 

0.2103 

1.11 

654.94 

Grid  search 

9 

0.2106 

1.11 

654.94 

Golden  section  Search 

9 

0.2102 

1.11 

654.94 

Hybrid 

12 

0.1737 

1.11 

654.94 

Losses 

9 

0.2102 

1.11 

654.94 

Loadability 

11 

0.4300 

1.39 

827.61 

maximization  (method  6)  gives  higher  losses  in  some  cases 
(12-bus  and  33-bus  test  systems). 

2.  Here  it  is  also  need  to  be  noted  that  the  capacitor  size  in  case  of 
loadability  maximization  (method  6)  have  also  been  signifi¬ 
cantly  increased  as  compared  to  other  methods.  However  in 
case  of  capacitor  placement  based  on  weakest  voltage  bus 
(method  5)  is  much  lesser  than  any  other  method. 


3.  The  voltage  profile  improvement  in  case  of  loadability  max¬ 
imization  (method  6)  was  found  better  than  the  other  methods 
in  all  test  cases,  as  shown  in  Fig.  2  (only  two  bus  test  system  are 
presented  in  the  result). 

4  To  visualize  the  impact  of  capacitor  placement  on  power  loss 
reduction,  voltage  profile  improvement  and  maximum  load- 
ability  improvement,  the  formulae  defined  in  Table  1  is 
calculated  and  plotted  in  Fig.  3. 

From  Fig.  3,  it  can  be  seen  that  the  losses  have  been  reduced  in 
case  of  methods  1-5,  however  the  method  6  can  increase  the 
system  losses.  However  in  terms  of  VPI  and  maximum  load- 
ability  improvement,  method  6  was  found  better  than  the  other 
methods. 

5  To  observe  the  maximum  amount  of  load  (kV  A)  the  system  can 
sustain,  the  load  on  each  bus  in  12-bus  test  system  in  presence 
of  capacitor  is  increased  till  the  allowable  voltage  limit  is 
reached  (0.95  <  Vbus  <  1.05)  and  the  results  are  summarized 
in  Table  4.  Here  it  can  be  seen  that  using  the  method  of 
loadability  maximization,  the  system  can  sustain  1.39  times  of 
base  load  (i.e.  1.39  x  595.4  kV  A=  827.606  kV  A),  however  other 
methods  can  sustain  only  1.1  times  of  base  load  (i.e.  1.1  x 
595.4  kV  A=  654.94  kV  A). 

6  One  of  the  important  factors  that  are  needed  to  be  considered 
in  capacitor  placement  is  reactive  line  current.  The  placement 
of  capacitor  may  increase  the  reactive  loading  on  some  of  the 
lines,  as  shown  in  Fig.  4.  It  is  interesting  to  note  that  the 
reactive  line  loading  on  some  of  the  lines  have  been  increased 
in  case  of  weakest  voltage  bus  approach  (method  5)  and 
maximization  of  system  loadability  (method  6)  as  compared 
to  base  case  and  minimization  of  loss  approach.  The  amount  of 
reactive  current  can  be  reduced  by  placement  of  multi¬ 
capacitor  units  in  the  system. 

From  the  above  discussion  it  can  be  concluded  that  the 
capacitor  placement  and  sizing  is  based  on  different  parameters 
including  minimization  of  power  losses,  voltage  profile 
improvement,  maximization  of  system  capacity  and  line  limit 
constraint.  Another  major  factor  that  is  needed  to  be  considered 
is  capital  cost,  the  capital  cost  will  include  capacitor  cost, 
protection  system  used  and  cost  of  constructing  new  lines  (if 
needed).  Methods  1-4  of  shunt  capacitor  placement  in  the 
system  can  well  handle  the  minimization  of  power  losses 
problem  considering  the  line  limit  constraints,  with  additional 
benefits  of  maximization  of  system  loading  and  voltage  profile 
improvement. 


5.  Conclusion 

This  paper  has  presented  a  very  detailed  overview  of 
optimum  shunt  capacitor  bank  (SCB)  placement  in  distribution 
system.  In  literature,  analytical,  numerical  programming,  heur¬ 
istic  and  artificial  intelligent  based  techniques  have  been 
proposed.  The  paper  has  also  presented  a  comparative  study 
of  six  different  techniques  of  optimum  capacitor  placement 
based  on  minimization  of  power  losses,  weakest  voltage  bus 
and  maximization  of  system  loadability.  The  results  show  that 
the  optimum  capacitor  placement  based  on  minimization  of 
power  losses  helps  in  reducing  the  reactive  current  component 
in  total  12R  losses,  in  addition  to  fulfill  the  line  current 
constraint.  The  other  advantages  including  voltage  profile 
improvement  and  maximization  of  system  loadability  are  also 
achieved  on  small  scale  using  power  loss  reduction  technique. 
Hybrid  method  (weakest  voltage  bus  and  minimization  of 
power  losses)  and  maximization  of  system  loadability  must  be 
utilized  carefully,  such  approaches  may  results  in  increasing  the 
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Fig.  4.  Comparison  of  line  current  due  to  loss  reduction  method,  hybrid  method  and  maximization  of  system  loadability.  (a)  12-Bus  Test  System,  (b)  30-Bus  Test  System, 
(c)  33-Bus  Test  System,  and  (d)  69-Bus  Test  System. 
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Table  A.1 

Proposed  Algorithm  for  Single-DG  Units  Placement 

Steps 

Disscussion 

Initialization 

■  Initialize  PSO  parameters. 

■  A  particle  i  among  initial  populations  (d)  is  a  one  or  two  dimension  position  vector  x„  representing 
shunt  capacitor  bank  position  (Qscb)  and  shunt  capacitor  sizes  (Qscb)  or  both. 

*i=  [Qscb  .  Qscb] 

Generate  Initial  Solutions 

Calculate  the  best  fitness  function  among  initial  solution  of  x,.  considering  the  constraint  given  in 

Eqns.  (3  and  4). 

Update  Initial  pbest  and  gbest 

For  Iteration =1,  find  initial  pbest  and  gbest • 

Generate  new  set  of  positions 

■  Update  velocity  and  position  vectors  using  Eqns.  (A-l )  and  (A-2). 

■  Check  the  feasibility  of  the  new  generated  particles  xi,  given  in  Eqns.  (3  and  4). 

Update  Initial  pbesc  and  gbest 

■  Calculate  new  fitness  function  of  updated  particles. 

■  Update  pbest  and  gbest. 

Stopping  Criteria 

Update  Iteration 

If 

Iteration  =  =  Iteration(maX) 

End 

reactive  line  current  in  some  of  the  lines  if  constraints  are  not 
defined  properly. 
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Appendix 


Particle  swarm  optimization  (PSO)  was  introduced  by  [131  ]  to 
solve  the  optimization  problem,  based  on  social-psychological 
metaphor  behavior.  A  particle  i  is  represented  as  x,-=(x,i,  x,2,  xi3, 

Xj4,  . ,  Xid  ).  The  position  associated  with  the  best  fitness 

pbesti=(pbestn,  pbesti2,  pbesti3,  pbestl4  . pbestid )  is  consid¬ 

ered  as  its  current  best  position.  Here  d  is  a  total  number  of  initial 
particles.  The  overall  best  position  of  the  population  associated 
with  the  current  overall  best  fitness  value  gbest  is  recorded.  The 
rate  of  the  position  of  ith  particle  is  represented  as  Vj=(vn,  va,  Va, 

12,-4, . vm).  During  the  iteration  procedure,  the  velocity  and  the 

position  of  ith  particles  are  updated  according  to  the  Eqs.  (A-l)  and 
(A-2): 

v,/+1)  =  wx  v,/>+Ci  xr,x  (pbestid-xid)+c2  x  r2  x  ( gbest-xid ) 

(A-l) 


xW(t+1)=xid(t)  +  vid(t+1) 


(A-2) 


where  t  is  number  of  iterations,  C]  and  c2  are  constants  (0.7),  r,  and 
r2  are  random  numbers,  w  is  inertia  weight  given  by  Eq.  (A-3). 


W  —  Wmax  +  ' 


iter, r 


-  x  t 


(A-3) 


where  wmax  and  wmin  is  0.9  and  0.4  respectively,  c2  and  c2  are 
constants  (ci  =  c2=0.7). 

Table  Al  Shows  the  complete  steps  of  the  proposed  algorithm 
for  optimum  single  unit  capacitor  placement  and  sizing. 

In  this  paper,  PSO  is  used  to  achieve  the  desired  fitness  function 
(minimization  of  power  losses  or  maximization  of  system  load- 
ability).  In  the  present  case  of  shunt  placement  and  sizing,  the  ith 
particle  (x,-)  is  a  two  dimension  vector  (Qscb.  Qscb).  representing 
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